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To clarify the conflict between the theoretical predictions and experimental measurements of 
the inclusive J/ifj production in T decay, We consider the a a order color-singlet contributions of 
processes T — > J ftp + gg and T — !> J ftp + gggg. Both the branching ratio and J ftp momentum 
spectrum are calculated, and the branching ratio (4.7 x 10 -4 ) is larger than the leading-order 
contribution (ai,T — > J ftp + ccg). Together with the QCD and QED leading-order contributions 
considered in our previous work, the color-singlet prediction of the branching ratio for the direct J/ip 
production is Br(T -> J/ip dilcct +X) = 0.901q;3? x 10 4 , which is still about 3.8 times less than the 
CLEO measurement. We also obtain a preliminary color-singlet prediction oiR cc = B ^(jr^Ij/^+xf 
and find the value 0.391q'20 ^ s much larger than the color-octet predictions, and suggest to measure 
this quality in future experimental analysis. 

PACS numbers: 12.38.Bx 13.25 Gv 14.40.Pq 



The existence of a hierarchy of energy scales: toq >> 
vtlqv >> vriQV 1 makes the heavy quarkonium system be 
an ideal laboratory to study both perturbative and non- 
perturbative aspects of QCD, where v, being assumed to 
be much smaller than 1, is the velocity of heavy quark in 
the rest frame of the heavy meson. And it is commonly 
believed that the nonrelativistic QCD (NRQCD) [l[ effec- 
tive theory provides a rigorous factorization formulism 
to separate the physics in different scales. The NRQCD 
not only covers the results of previous color-singlet(CS) 
model predictions, where at short distance the QQ can 
only be in the color-singlet configuration with the same 
quantum numbers as the corresponding heavy quarko- 
nium states, it also includes the contribution of QQ in 
the color-octet(CO) configuration at short distance. 

Despite of the impressive success of the NRQCD, the 
role of the CO is not well established yet, particularly 
in J/tp production case. The substantial theoretical 
progress in the next-to-leading-order(NLO) calculations 
shows that there is not a convincing mechanism to ex- 
plain the J/ip production data in varies experiments self- 
consistently yet. For J/^+cc+xJU and J/?A+A" non _ cE 
production in e + e~ annihilation at B-factories, the CS 
processes themselves can account for the cross sections 
when the NLO QCD corrections an d relativistic cor- 
rections [1, H[ are taken into account, which leaves very 
little room for the CO contribution [10|. For J/ip pro- 
duction from Z decay, recent NLO QCD correction cal- 
culation in the CS ll| gives just one-half of the experi- 
mental measurement and the CO contribution might be 
able to explain the other-half. The transverse momen- 
tum p t distribution of J/ip photoproduction and polar- 



ization parameters at HERA can not be well described 
by the CSjH, [H at QCD NLO, it seems that p t dis- 
tribution of J/ip photoproduction can be explaine d by 
the CO and CS contribution together at QCD N Loll5l: 
For J/ip hadroproduction, together with the CS 16l4l8l| 
contribution and the CO contribution [l9|, [2(| at NLO 
in a s , we can not describe the Tevatron results about 
the pt distribution of J/ip production and polarization 
simultaneously yet. 

In order to clarify such a puzzling theoretical situation, 
it is worth to further investigate some other J/ip produc- 
tion processes, one of which is the inclusive J/ip pro- 
duction in T decay. From the theoretical point of view, 
because T predominately deca ys i nto three giuons via bb 
annihilation, it is proposed [2lil22j| that in the rich-gluon 
final state environment abundant J/ip can be produced 
through cc pair in CO 3 Si configuration. Hence, the in- 
clusive J/ip production in T decay will be an another 
good probe to discriminate the CS and CO mechanism. 
And the present CO predictions of the branching ratio 
is B(T -> J/ip + X) = 6.2 x 1CT 4 (HI, S3 with about 
10% feed-down contribution from ip(2S) and 10% from 
XcJ [HI- While the correct CS result, which was overes- 
timated by about an order in magnitude before [ID, [2(| , 
is only 4.2 x 10 _5 [27]]. On the experimental side, the 
branching ratio for T — > J/ip + X has been measured 
by a few collaborations about twenty years ago [281 - 4301 ] . 
and recently, the more precise measurement carried by 
the CLEO Collaboration gave 31 



B(T -> J/i/j + X) = (6.4 ± 0.4 ± 0.6) x 10" 



(1) 



* Present address 



It can be seen that the CLEO result is in good agree- 
ment with the CO prediction, but the J/ip momentum 
distribution measured by CLEO [3l[ is much softer than 
the CO predictions [2l|, [22j. In a very recent work[23|. 
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it is found that the momentum spectrum can be signif- 
icantly softened after combining the NRQCD and Soft 
Collinear Effective Theory(SECT) in the kinematic end- 
point region. However, it yields a much smaller branching 
ratio. This may indicate that the CO processes do not 
contribute dominantly. 

In the CS case, the results at QCD and QED leading- 
order(LO)[13| is about an order of magnitude smaller 
than the experimental result. To investigate the CS 
contributions for the J/ip production in T decay more 
precisely, we consider the 0{oP s ) contributions through 
T — > J/ip + gg and T — > J/ip + gggg processes. It 
was refereed in Ref. [24j that this two processes had been 
crudely estimated and the branching ratio is a few of 
x 1CP 4 , which is comparable to CLEO data. Since the 
former calculation refereed is rough without any publi- 
cation and detail, it is necessary to give an exact result 
and a complete analysis of these two processes. To per- 
form the calculations, we employ the Feynman Diagram 
Calculation (FDC) package j32j j . 

According to NRQCD factorization approach, at the 
leading order of v\ and v\, the CS contribution to T ^ 
J/ip + X is expressed as: 

fff(T -> J/ip + X) = df(bb[ 3 S! , 1] -> 

C c[ 3 Sx,l] + X) x (Tp^S^iOfeS,)), (2) 

where dT is partonic partial decay width which can be 
calculated perturbatively. By dimension analysis, it is 
easy to derive out that the general expression of the par- 
tial width is written as 
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(3) 



culation, each diagram is summed up with it's possible 
partner diagram which is obtained by reversing the direc- 
tion of b-quark or c-quark line in fermion loop. Then the 
diagrams are divided into ten groups, and the representa- 
tive ones are shown in Fig.fl] and others can be obtained 
by exchanging the positions of the two final-state gluons. 
For diagrams in the same group, the amplitude equals 
to each other when ignoring color factor, thus only the 
d abc d abc piece in color factor will survive after all the di- 
agrams being summed up in each group. It is found that 
there are infrared divergences in the amplitudes of di- 
agrams in Al, A2, A4, A5 groups, at least in Feynman 
gauge. And the divergence terms in A1(A2) group cancel 
those in A4(A5), then the total amplitude is finite. The 
amplitude of each diagram in group A3 is finite individu- 
ally because in such diagrams there are no such a virtual 
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FIG. 1: The typical Feynman diagrams of five of ten groups 
in the CS T — > J/ip + gg process, and the other five typical 
ones can be obtained by exchanging the positions of the two 
final-state gluons. 



where r = m c /m,b is a dimensionless parameter and / is 
a process dependent function of r. The (T|C 1 ( 3 S' 1 )|T) 
and (Of ( 3 Si)) in Eq.[5] arc the nonperturbative matrix 
elements, which will be determined phenomenologically. 
To be consistent with our former work (27j . we keep the 
factor 3 ( 2 jv p explicitly. 

The 0(a®) T — > J/ip + gg process is of 36 one-loop 
Feynman Diagrams at leading-order. To simplify the cal- 



gluon which joints two on shell (anti) quarks or one on 
shell quark and one on shell antiquark. 

Before showing the result, we would like to address 
some non-trivial techniques treatment in the calculation. 
By applying the FDC package, the general expression of 
the Feynman amplitude for T(pi,ei) — > J/ip(j>2,£2) + 
ff(P3i £3) + 9{pa, £4) process is generated as: 



M = ei • e 2 e 3 • e 4 c 4 i + ei • £ 3 £ 2 • e 4 c 42 + ei • £462 ■ £ 3 c 43 + Pi • e 2 Pi ■ £3£i ■ £4Ci 9 + p 4 ■ £i(P4 • £2£3 • £4037 

+Pi ■ £3£2 • £4239) + P3 ■ ti[Pi ■ £l(£2 ' £3C21 +P4 ' £2^4 • £3Cq) + Pi ' £2£l ' £3C29 + Pi ■ £3£l ' £2C 3 l] 

+P3 ■ tl[Pi ■ £3£l • £4C 2 +Pi ■ £l£3 • £4C35 + P3 ' £4<£l • £3C27 + Pi ■ tlPi ■ £3Clo)] + Pz ■ £l[P4 ■ £2£3 • £4C 3 8 

+Pi ■ £3£2 • £4C40 + P3 ■ £4(£2 ' £3C22 + Pi ' £2P4 ' £3Cll) + P3 ' ^2^3 ' £4C 3 6 + P3 ' £4P4 ' £3212)] 

+P2 ■ £3(^4 ' £2£l ' £4Ci 7 +Pi ■ £l£2 • £4C 33 +P2 • £4^1 ' £2C 2 6 + P3 ' £2^4 • £lC 6 + P3 ' £l(p 3 ' £2C 8 (4) 
+Pi ■ £2C 7 ) + Pi ■ eiPi ■ e 2 c 5 ] + p 3 ■ e 4 (ei • e 2 c 2 5 + Pi • £iP4 ■ eaCi) + p 3 ■ e 2 (ei ■ £4Cis + p 3 ■ tiPi ■ £ic 2 ) 

+P3 ■ £l(£2 ' £4C 3 4 +P3 ■ tiPi ■ £2C 3 + P3 • £2j» 3 ■ £404)} + P2 • £4[P3 ' £2(fl • £3C28 + Pi ' £lP4 • £3014) 
+P3 ■ £l(£2 ' £3C24 +P3 ■ t2Pi • £3Cl6 + Pi ' £2P4 • £3215) + Pi ' £l(£2 • £3C23 + Pi ' £2P4 ' £3Cl3) 
+Pi ■ £2£l • £3C30 + Pi ■ £3£l • £2C 32 ] 

I 

where Cj, i — 1, . . . , 43 are the coefficients of the Lorentz structure and all the loop diagrams contribute to their 
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values. From the general expression in Eq.Q, it is easy 
to understand why the tensor reduction procedures are 
very complicated and will generate complicated results 
for the coefficients Cj. 

These complicated results of Cj may contain fake pole 
structures and cause big number cancellation problem, 
finally will spoil numerical calculation in limited preci- 
sion case. In fact, it really took place in our numerical 
calculation. We found it is impossible to control the can- 
cellation of the big numbers and to obtain correct results 
in double precision FORTRAN calculation. In quadruple 
precision FORTRAN calculation, we can obtain the cor- 
rect results by introducing cut-conditions to control these 
fake poles in phase space integration. To demonstrate 
that our treatment is suitable, we define a cut condition 
parameter in phase space integration as 

_ T(|M| 2 = 1, with cuts) 

~ T(\M\ 2 = 1, without cuts) ^ ' 

and calculate the partial decay width numerically with 
different cut condition parameter. The results are: 

T = a(1.440 ± 0.001) x 10~ 10 for A = 1 - 0.04 
T = o(1.446 ± 0.001) x 10" 10 for A = 1 - 0.005 
T = a(1.450± 0.001) x 10" 10 for A = 1 - 0.0005, 

(6) 

where all the calculations are under control in quadruple 
precision FORTRAN calculation, and a is just a con- 
stant number. We have tried to do the calculation with 
A = 1 — 0.00005 and found a divergent result. It means 
that the big number cancellation will lose control even in 
quadruple precision FORTRAN when A is too close to 1. 
Therefore, our following calculations for T — > J/ip + gg 
are based on the cut condition parameter A = 1 — 0.0005. 
Moreover, to check the gauge invariance, in the expres- 
sion of Eq. (j4]) we replace the gluon polarization vector 
(or £4) by its 4-momentum p^ (or p^) in the final numer- 
ical calculation. Definitely the result must be zero and 
our results reproduce it. 

Now, we proceed to present our results. Since the two 
interior gluons can be on shell simultaneously in this pro- 
cess, the amplitude will be complex- valued. And we use 
the superscript 'Im" to denote the contribution of the 
real process T — > 3g followed by gg — > J/ip + g. Setting 
m c = m J/4 ,/2 = 1.548 GeV and m b = m T /2 = 4.73GeV, 
which corresponds to r = 1.548/4.73 = 0.327, we get 
f gg (r) = 2.07 and f^(r) = 0.741 which is about 1/3 of 
the total. To show the dependence of / on r, we also list 
some of the numerical results of f(r) in Tab. I, where r is 
in the range of 0.275 < r < 0.381, which is obtained by 
fixing the value of rrib and varying m c from 1.3GeV to 
1.8GeV GeV[25j|. For comparison we also list the results 
of fecgir) for the 0{a 5 s ) T — > J/ip + cc + g process. It can 
be seen that both the value of f gg {r) and that of f^(r) 
do not change sharply, when r goes from 0.275 to 0.381, 
and this behavior is quite different from what happens to 
fccg{r)- Also the ratio of fgf{r) to f gg (r) changes very 
little with r. 



TABLE I: The values of f(r) for J/ip + cc + g (f C c g ), J/ip + gg 
(f gg and fg^),J/ip + gggg {fig) production in T decay with 
different inputs of r = — . 



r 


J ccg \ 1 j 


faa(r) 
J 99 v / 


J 99 V / 




0.275 


0.904 


2.94 


1.02 


1.35 x 10" 2 


0.296 


0.567 


2.54 


0.892 


1.08 x 10~ 2 


0.317 


0.345 


2.21 


0.786 


0.880 x 10~ 2 


0.327 


0.269 


2.07 


0.741 


0.800 x 10~ 2 


0.338 


0.202 


1.94 


0.696 


0.721 x 10~ 2 


0.361 


0.105 


1.68 


0.612 


0.585 x 10~ 2 


0.381 


0.055 


1.49 


0.547 


0.490 x 10~ 2 



There are 216 Feynman diagrams in the CS T — » 
J/ip+gggg process, and the typical one is shown in Fig. 2. 
It is a tree process without infrared divergence, and the 
numerical results are calculated straightforwardly with 
the help of FDC package. When r = 0.327, we get 
fi g { r ) = 0.8 x 10~ 2 , which is more than two orders less 
than / ffS (0.327). Some other numerical results of f± g (r) 
for 0.275 < r < 0.381 are also listed in Tab. I. Like f gg (r), 
the function fi g (r) also does not dependent on r seriously, 
but its value is too small comparing to the values of the 
/ functions of the other processes. One possible reason 
is that the five-body phase space is much smaller than 
the three-body phase space. 

Besides f gg (r)(fi g (r)), the partial decay width T(T — > 
J/ip + gg{^g)) also dependents on the choice of the val- 
ues of the two NRQCD long-distance matrix elements, 
the coupling constant a s and the b-quark mass mj. The 
value of (Cf ( 3 S"i)) ^ 3{J/i)\0( 3 Si)\ J/tp) can be ex- 
tracted from J/ip decay into e + e~ by using the upto a s 
order result 

r(j/v> e+ e -) = 2 -^(i - ^widfSxM. 

6m* An 

(7) 

Using a = 1/128, m c = 1.548GeV, a s {2m c ) = 0.26, 
T(J/ip -> e+e") = 5.54keV@, we get (0f( 3 5i)) = 
1.25GeV 3 . And (T|C( 3 S' 1 )|T) = 2.92GeV is be deter- 
mined in a similar way with rrib = 4.73GeV, a s (2rrib) = 
0.18, T(T -> e+e') = 1.29keV0. The uncertainty from 
the choice of the renormalization scale is quite large since 
there are two typical energy scales rrib and m c in the cal- 
culation. By choosing fi — 2m c , we find 

rs = T(T -> J 14, + gg) + T(T -> J/iP + Ag) 

= 9.1xlO~ 3 keV (8) 

It corresponds to 

B 9 = B(T -> J/ip + gg) + B(T -> J/ip + 4g) 

= 1.7xl0- 4 (9) 

which is coincident with the rough result mentioned in 
Ref.Jlij]. However, the branching ratio becomes much 
smaller and is only 2.32 x 10~ 5 when choosing /i = 2rn,b- 
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FIG. 2: One of the 216 Feynman diagrams for the CS T 
J/i> + 9999 process. 



This is because the processes are at order. Note that, 
seemingly, the branching ratio also strongly dependent 
on mb as m^ 5 , but in fact the dependence is m^ 3 be- 
cause of the dependence of the nonperturbative matrix 
element (T|C( 3 5i)|T) on nib from it's phenomenological 
determination. To obtain the above numerical results for 
the branching ratio, the experimental measurement on 
the T total decay width Ty = 53keV 33] is used. 

Since the J/tp is produced through the three-gluon 
decay channel of T in these two processes, It is natu- 
ral to normalize the partial width to the decay width of 
T — > ggg, which at LO in a s is given by 

T(T ggg) = ^j^ 2 - 9)(T\0( 3 S^T) . (10) 

Then branching ratio is expressed in an alternate form 



B° = r« or x B(T ggg), 



where 



1 Nor 



81(/ flfl (r)+/ 4fl (0)ag(Or( 3 Si)) 
20{2N c ) 2 m 3 b (7T 2 -9) 



(11) 



(12) 



and B(Y —> ggg) = 84% is obtained by assuming B(Y — > 
ggg) « B(T -> light hadron (LH)) a 33]. To calculate 
the branching ratio in this way is equivalent to determine 
a 3 (T\0( 3 Si)\T) from LH decay of T and can reduce the 
uncertainties from a s . And our following results are all 
calculated based on Eq. [11] . 

The numerical results of "f(r)" for each decay process 
are presented in Tab.l and it show that the values of 
f( gg )(r) changes slowly when r goes from 0.275 to 0.381. 
The Feynman diagrams in Fig.fl] indicate the process 
T — » J /if) + gg can be viewed as T — > gg(*'g(*' b followed 
by 3 ( *- ) 3 ( *' ) -> J/ip + g- Then the normalized F^ or can 
cancel part of the contribution at m& scale, so similar to 
what is done in Ref. [2l| we choose the scale of a s to be 
2m c . The theoretical uncertainties of T — » J/ip + gggg 
can be analyzed in the same way. 

By setting the default parameter choice: nib = 
4.73GeV, r = 0.327, (Of( 3 S'i)) = 1.25GeV 3 and 



a The contribution of T — > 7* — > qq is excluded. 

b means the gluon can either be virtual or real. 



J/v+gg, 




FIG. 3: The normalized partial width of the 0(a 6 s ) T -4 
J/xp + gg process as function of J/ip momentum pj/^. The 
solid line is the total result, and the dashed line is the contri- 
bution from the imaginary part in the Feynman amplitude. 



a s (2m c ) — 0.26, we obtain 
B 9 



r« or x B(T ^ ggg) = 0.47 x 10" 



(13) 



Using the same inputs to re-estimated the result in 
Ref. 27] and adding up it with the contribution in Eq.[l3 . 
we obtain the total CS singlet prediction 



B(T -4 J/ip + X) = 0.90 x 10" 



(14) 



where the total contribution from the 0(a 6 s ) J/ip+gg and 
J/ip + MR 9 processes is as important as those calculated 
in Ref. [27]. It is clear that the uncertainties are from 
the b quark mass mj, the scaleless functions f gg (r) and 
f4g(r), and the choice of the scale of a 8 . To estimate 
the uncertainty, we used nib = A.&GeV, r — 0.296 and 
/i = 2m c for upper bound; mb — 4.9GeV, r = 0.361 and 
fx = 2m c for lower bound, then the branching ratio is 
represented as: 



B(T J/ip + X) 



0.90ig-j? x 10" 



(15) 



Furthermore the total branching ratio turns to be a 
much smaller value 6.3 x 10~ 5 (5.2 x 10~ 5 ) by choos- 
ing the scale to be 2nib (2^/7n&m c ) and a s (2nib) = 0.18 
(a s (2y/m c mb) — 0.21) and keeping the other parameters 
the same as for the central value of the branching ratio. 

The experimental result in Eq.([T|) includes the feed- 
down contributions of XcJi which are < 8.2,11,10 per- 
cents for J = 0, 1, 2 respectively, and 24% feed-down con- 
tribution of ip(2S). Removing the feed-down contribu- 
tions, the branching ratio of direct J/ip production in T 
decay would be 



B(T -4 J/V-dhcct + X) = 3.52 x 10" 



(16) 
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FIG. 4: The normalized partial width for the CS J/tp pro- 
duction in T decay as function of J/tp momentum Pj/j,. The 
solid line is the total result, the dashed line is the contribu- 
tion of the 0(t*s) T — > J/ip + gg process. The dashed-dot line 
is the lOOx the contribution 0(a & a ) T — > JJtb + gggg. The dot 
line is the contribution calculated in Ref. 27], which includes 
the 0(a 5 s ) T J/tp + ccg and 0(o?a 2 s ) T -> J/ip + gg and 
T — > J/ip + cc processes. 



Which is about 3.8 times larger than the current CS pre- 
diction. 

For the J/tp momentum spectrum, the normalized de- 
cay widths denned in Eq.([i"2"|) are used to present the re- 
sults with default parameter choice. Both the total result 
(solid line) and the imaginary part contribution (dashed 
line) for 0(a 6 s ) J/tp + gg process are shown in Fig. [3] 
with very similar shape. A summarized CS contribution 
to the pj/jj, distribution of the normalized decay width 
is shown in Fig. [4] . And we find the peak of total result 
curve is at pj/$ — 2.7GeV, which is a little larger than 
that of the CLEO measurement [3l| . 

It is found that the J/tp production in associa- 
tion with cc pair is an important mechanism for J/tp 
electroproduction0,|3| in the Belle experiment. And the- 
oretically, the contribution of pp — > J/tp + cc+X process 
to J/ijj hadroproduction at the Tevatron is also found to 
be non-ignorable [35i l3a| . The ratio of J/tp production in 
association with cc pair to J/tp plus anything may also 
be a good probe to reveal the J/tp production mecha- 
nism in T decay and to clarify the conflict between the 
CLEO measurement and theoretical prediction. Choos- 
ing a s (2m c ) = 0.259, we give the CS prediction for the 
ratio R rr 



R c 



B{T -» J/tP + cc + X) 
B(T ^ J/tP + X) 



0.39: 



0.20' 



(17) 



where the center, upper and lower bound values corre- 
spond to r = 0.327, 0.296 and 0.361 respectively, and the 
associated charmed particles process includes the 0(a^) 
T — > J/tp + ccg sub-process, which is dominant, and 



FIG. 5: The normalized partial width of T decay into J/tp + 
cc + X process as function of J/tp momentum pj/^ (dashed 
line) and that for J/tp + X production (solid line). 



0{a 2 a 2 s ) T — > 7* —> J/tp + cc sub-process. On the con- 
trary, the CO prediction of R cc is only at the level of 
1%|21|, which is quite different with the CS prediction. 
Unlike the branching ratio, the theoretical prediction of 
R cc only depends on r and a s , which results in a rela- 
tively small uncertainty. Particularly, if we drop the con- 
tribution of QED part, R cc is just proportional to a s . In 
Ref. [1, [j| , the authors find the enhancement of the NLO 
QCD corrections is large in e + e~ — > J/tp + cc process. It 
indicates that the result in Eq. [T7] is only a very prelimi- 
nary result and to get a more solid predictions the contri- 
bution of the NLO QCD corrections to T — > J/tp + cc + g 
process should be taken into account. Calculating the 
NLO QCD corrections to T — > J/tp + cc + g process is 
beyond the scope of this work and will not be discussed 
here. In the end, the J/tp momentum spectra for the as- 
sociated process and non — cc process are given in Fig. [5] 
for comparison. 

In summary, in this work, we calculate the 0(a 6 s ) CS 
contribution of T — > J/tp + gg and T — > J/tp + gggg pro- 
cesses to the inclusive J/tp production in T decay. The 
branching ratio is estimated in two ways. In the first 
way, the numerical results of partial width and branch- 
ing ratios are all evaluated directly. And we find its re- 
sult is more close to the experimental data when the 
values of the parameters are properly chosen, which is 
also coincident with a rough estimated result mentioned 
in Ref. [24]. However, the uncertainty of this way is 
very large and the branching can be in a wide range of 
2.3 x 10~ 5 ~ 1.7 x 10~ 4 . In the second way, the branching 
ratio is calculated by using the normalized decay width, 
which seems more reliable. After combining the present 
result with the contribution calculated in our previous 
work[27j|. we find now the total CS prediction is about 
0.92 x 10~ 4 , which is still about 3.8 times less than the 
experimental value 3.2 x 10~ 4 for direct J/tp production 
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given by the CLEO Collaboration [H|. We also calculate 
the J/ip momentum spectrum and find the peak of the 
color-singlet curves is close to that of the CLEO result, 
although being a little larger. Besides the the branch- 
ing ratio and the J/ip spectrum, we also study the ratio 
R cc = B(T -> J/ip + cc + X)/B(T -> J/ip + X), and 
find the CS prediction is much larger than that of the 
CO. Since the associated charmed meson process can be 
measured separately, a analysis of the ratio is expected to 
perform in the CLEO, Babar or Belle experiments. Now 
there is still large discrepancy between the CLEO results 
and NRQCD predictions. There are two points to be ad- 
dressed: first the J/ip production mechanism is not well 
understood yet, and the existence of the CO mechanism 
is still under debate; second the higher QCD corrections 
are not included completely. Therefore, to understand 



the J/ip production mechanism in T decay and more- 
over in pp collisions at the Tevatron. Further theoretical 
and experimental work are necessary. 
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